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to define a non-Gaussian density operator. 
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The Duffing oscillator is a typical anharmonic oscillator whose classical aspect is rel- 
atively well- understood |jl||. It has been studied in many fields of physics from classical 
mechanics to quantum mechanics, especially as a toy model for one-dimensional quantum 
field theory. In particular, more physical interests have been focused on the application of 
various perturbation techniques both as a classical and as a quantum model |2J. 

In this paper we study the quantum Duffing oscillator, a time-independent quantum 
anharmonic oscillator, 
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Even though one can find in principle the exact quantum states by solving the time- 
independent Schrodinger equation, the results of most of attempts are minimal. 

In this Letters to Editor we shall propose a mean-field approach to the quantum Duffing 
oscillator, which has the following points different from other approaches: It takes into 
account a physically meaningful frequency that is quite close to the exact classical one. 
There appear no secular terms at least at the order of mX for the coupling constant up 
to a critical value. A kind of renormalization can be used at the order of (mA) 2 . The 
motivation for this approach is based on two observations: both the extremization of the 
Dirac action and the quantum Liouville equation at the quadratic order with respect to the 
vacuum state of a Fock space lead to the same mean-field equation. We find a complex 
solution whose frequency is the same as the exact classical one at the order of mX, construct 
the operators on the Fock space that generate the exact quantum states perturbatively, and 
define a non-Gaussian density operator. 

A general Fock space is constructed by the basis (h = 1) 
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* = up + vq,a = u*p + v*q (2) 



where (u*v — uv*) = i follows from the usual commutation relation 



a, a* 



1. u and v 



depend on time explicitly, as will shown later. We extremize the Dirac action 



T I dt(*,t\i^-H\*,t). (3) 
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The variation of Eq. (||) with respect to \^?,t) leads to the time- dependent Schrodinger 



equation, i-^ \ty , t) = H \^ , t) . The exact quantum state of the form 



(4) 



gives the null action. As an approximate quantum state we use the time-dependent number 
state a^a \n,t) = n \n,t) of the Fock space. Let v = —mil, n* = mil, and n = mu* . Then, 
the mean energy of the ground state 
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(0, t\ H\Q,t) = — + —u*u + —{u*u) 2 , 
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acts as a Hamiltonian. The Hamiltonian equations is equivalent to a second order mean-field 
equation 



u + lo 2 u + 3X(u*u)u = 0. 
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We may derive Eq. (|j) from a different argument, the quantum Liouville equation 
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We expand the Hamiltonian operator in terms of the creation and annihilation operators 



H = H 2 + m\H 4 + (0, t\ H |0, t) 



(8) 



where H 2 and H± are quadratic and quartic in a* and a. Then, a^ and a satisfy the approx- 
imate quantum Liouville equation 
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A complex solution to Eq. (0) is found 
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where fl satisfies a cubic equation 
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There is at least one real positive root 
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It is remarkable that the mean field equation takes already the same frequency at the order 
of m\ as the exact classical frequency Q, which is different from those in Refs. 0Q. 

We construct the operators of the quantum Liouville equation (0) in a perturbative way 



it = a t + ^( mA ) n 4 +ll i = a + ^(mA)"B 2 n + i 
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where 
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For a weak coupling constant we solve Eq. (|7]) order by order in mX. The lowest order 
equation is Eq. (0), which has already been solved. At the next order, we obtain 
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The first bracket on the right hand side is satisfied provided that a' and a satisfy Eq. @ 
and all the derivatives act only on the coefficient functions b 3 . but not on a) and a. The 
second and third brackets lead to the inhomogeneous equation 
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For a solution of the form 



b 3 := (b 3 ,b 3 ,blb 3 ) 
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the homogeneous equation becomes a system of differential equations with constant coeffi- 
cients 
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where a = 4n ^ 3 . All the eigenvalues of Eq. (|18|) are real constants up to the critical coupling 



constant a r 
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0.1365. For the coupling constant larger than A c , a secular term of an 



exponential function begins to appear. This implies that for the strong coupling limit the 
Fock space basis may not be suitable for constructing bounded operators such as A* and A 
that generate the quantum states. On the other hand, for the weak coupling limit (a << 1), 



numericals result show that b 3 ~ (c^e 
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). For a = 0, the 



result is exact with the eigevalues (— 4f2, —20,0,2^). The freedom in choosing the norms 
of eigenvectors of Eq. (|l^) may be used to keep the same equation as Eq. @, a kind of 
renormalization, even at the order of (mX) 2 . Thus, the operators, Eq. (|i~3"D, are bounded 
operators whose coefficients are proportional to A 



in 2 ' 



Finally, we discuss on the density operator defined as 
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-ciqA^A 



(19) 



The density operator, Eq. (|I~9"D, satisfies manifestly the quantum Liouville equation (0) by its 
construction at the order of mX. It would be worthy to compare it with the Gaussian-type 
density operator that was constructed for a time-dependent Duffing oscillator ||, because 
the operators A^ and A can also work for the time-dependent Duffing oscillator (in fact, the 
complex solution (^) is an adiabatic one in this case) and, furthermore, Eq. ([1~9|) already 
went beyond the quadratic order for the Gaussian-type. 

The details of this approach and some problems such as the divergence structure of the 
operators A^ and A and the relation with the renormalization group approach |||| will be 
addressed in a future publication fTj. 
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